The (In) stability of the Bianchi IX Dynamics: 
The Geodesic Deviation Equation in the Finsler Spaces 



M. Di Bari 1 '* and P. Cipriani 2 ' 1 
1 Dipartimento di Fisica, Universita di Parma, viale delle Scienze, - 43100 - Parma. 
2 Dipartimento di Fisica, Universita "Tor Vergata", via della Ricerca Scientifica, 1 - 00133 - Roma. 

(June 16, 1998) 



00 

On 
OS 



3 

OS 



> ■ 

<n : 

(N . 

o ■ 
■ 

o : 

00 . 
Os ■ 

"o ! 

i < 
5h . 

bo 



We explore the dynamical stability of the minisuperspace 
Hamiltonian of the Bianchi IX cosmological models, giving 
a gauge-invariant and unapproximated description of the full 
continuous dynamics, achieved through a geometrical descrip- 
tion of the equations of motion in the framework of the theory 
of Finsler Spaces. The numerical integrations of the geodesies 
and geodesic deviation equations show clearly the absence 
of any "traditional" signature of Chaos, while suggesting a 
chaotic scattering dynamics scenario. 
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The characterization of Chaos in the framework of 
General Relativity is still an open question, receiving re- 
cently an increasing attention [Q . The issue is related to 
the well known ambiguities hidden in any gauge theory, 
and in particular to the freedom in the choice of the co- 
ordinate system. A gauge invariant description of Chaos 
must be invariant under any such choice, in particular 
with respect to rescaling of the "time" parameter. The 
problem is mainly related to the choice of an invariant in- 
ternal time which must be adopted in describing the 
evolution of a general relativistic dynamical system. The 
applications to the Bianchi IX dynamics of the standard 
approaches to Chaos borrowed from classical theory of 
dynamical systems have, in general, revealed their inad- 
equacy to satisfy this requirement. This has been out- 
lined in the past |j| , in particular for what concerns the 
dependence of the Lyapunov characteristic numbers on 
the choice of the time variable. 

The main works devoted to the study of Bianchi IX dy- 
namical behaviour are essentially of two kinds. On one 
hand the first studies of Belinski et al. |i) and the subse- 
quent papers || , till the last of Cornish and Levin || are 
based on a reduction of the continuous flow to discrete 
maps (i.e. Gauss and Farey maps) which mark the more 
relevant events in the evolution. In this case the time 
is a discrete variable which increases once a particular 
event occurs, i.e. when there is a bounce on the poten- 
tial walls in the minisuperspace framework |F3] or, in the 
other picture, a crossing to a new Kasner epoch and/or 
era Q . These events provide a natural clock for the evo- 
lution of the system. The discrete approximation of dy- 
namics by maps is usually adopted also in classical me- 
chanics [p| , where the definition of time is not questioned 
and the interval between two successive increments in 



the maps' variables, is roughly constant. This is not 
the case for Bianchi IX models, as the lapse of time de- 
pends on the definition of time itself, unless an invariant 
choice is adopted. In B, with the aim of giving an intrin- 
sic description and relying on the observer-independent 
properties of fractals in phase space, the authors argue 
a chaotic behaviour of the full dynamics by calculating 
the strange repellor's Lyapunov dimension, the topologi- 
cal entropy and the multifractal dimensions of the Farey 
map. 

However, even if the maps reproduce the main character- 
istic of the flow in the Kasner regime, they constitute only 
a discrete approximation of the actual field equations for 
the Bianchi IX models; so, it is reasonable that some 
efforts have been devoted to analyse the continuos dy- 
namical flow. Within the framework of a gauge-invariant 
description of motions and of their possibly chaotic prop- 
erties, the geometrical description of the dynamics offers 
a natural tool, as it gives a metric in terms of which also 
an invariant time variable is defined. The first attempt 
towards this kind of approach dates back to Chitre M, 
which reduced the motion to a geodesic flow on a hyper- 
bolic manifold, in order to advocate the mixing proper- 
ties of geodesies on a compact Riemannian manifold with 
constant negative curvature [flOf| . However, some approx- 
imations need to be assumed, which are ultimately equiv- 
alent to those involved for the maps. In the hope of avoid- 
ing any approximation, some authors [pd] adopted the 
Jacobi metric which derives from the Maupertuis princi- 
ple and which has been exploited in classical mechanics, 
in order to understand the sources of Chaos in dynamical 
systems jl2|. The internal time is given in this case by 



dsj =V2[E- U{x 1 )] dt 



(1) 



where E is the total energy and U(x l ) is the potential 
energy depending on the coordinates {x 1 } alone. The 
manifold is again Riemannian but the curvature is not 
constant and negative. The features of the motion are 
now ascribed to the average behaviour of the Ricci curva- 
ture. A clear detection of a predominantly negative Ricci 
curvature, which goes to — oo whenever there is a bounce, 
would indicate the presence of instability, at least in a lo- 
cal sense rather than as a global property of the motion 
p3| . Neverthless, in a general relativistic context, it is 
E = 0, because of the covariant character of the theory, 
so T = -U and dsj = V^Tdt = -V2Udt. Therefore the 
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kinetic term is not positive definite and, even imposing 
by hand that dsj = y/2\T\dt = V2\U\dt, the relation be- 
tween the time t and the "Jacobi parameter" sj becomes 
meaningless when T = — U = 0, which leads to singular- 
ities in the sectional, Ricci and scalar curvatures. So, 
although a Jacobi metric can be introduced without any 
approximation, in order to analyse the curvature prop- 
erties, it is necessary to cut the trajectory when these 
singularities occur. For this reason, the Jacobi metric 
cannot be properly used |l4j . As noticed in p5[ , here the 
zero kinetic energy surface no longer constitutes the ex- 
ternal boundary of the region allowed to motion in config- 
uration space, so it can be crossed infinitely many times 
by the trajectory, with different incident angles. More- 
over, as Sj runs essentially only near the potential walls, 
this description mimics again the reduction of the flow 
to a discrete map |p|Jl5l|. In such a way, both the maps 
and the Riemannian pictures, base their results about 
stability on substantially equivalent approximations. In 
conclusion, an unapproximated method, which can give a 
full description of the eventual instability of the Bianchi 
IX dynamics, is still lacking. 

Aim of this work is to provide such a description in the 
framework of the geometrodynamical approach to Chaos, 
generalizing the Jacobi description of motion in order to 
avoid any singular behaviour (see Ji6|]). We explicitly 
analyze the stability properties by numerically integrat- 
ing the geodesies and geodesic deviation equations in this 
more general geometrical setting. 

The generalization of Riemannian manifolds is accom- 
plished using the theory of Finsler Spaces jlTj , in which 
the line element, i.e. the invariant time parameter, and 
all the geometrical quantities (in particular the metric 
and the curvatures) depend on both coordinates and ve- 
locities. The Finsler geometrodynamical approach to 
Chaos has been already exploited in other dynamical 
systems. In we compared the tangent dynamics 

with the geodesic deviation equations of the Jacobi and 
Finsler metrics to highlight differencies and analogies, 
taking as paradigmatic classical dynamical system, the 
Henon-Heiles Hamiltonian. This system has been widely 
investigated in Pg| , where we find a strong evidence of 
correlation between geometric properties of the Finsler 
manifold and regular or chaotic behaviour of single or- 
bits. In the Finsler metric, the natural definition of the 
invariant interval of time is given by 



ds = ds F = Ciq^q^dt (i = l,...,N), 



(2) 



where C(q l , q 1 ) is the Lagrangian of the system and N is 
the number of degrees of freedom of the system. In ]l9| ] 
we developed the mathematical formalism to analyze the 
instability properties of very general Lagrangian systems 
with N degrees of freedom and performed some prelimi- 
nary steps towards a gauge- invariant description of insta- 
bility for the Bianchi IX dynamics. In order to consider 
the definition (0) as a good one and to derive from it a 



metric tensor, the Lagrangian must fulfill, among others, 
the condition to be sign definite (e.g., positive). This can 
be obtained quite simply, remembering that it is always 
possible to add to the Lagrangian the total derivative of 
a function of coordinates, without changing the dynam- 
ics. As the condition C > can be fulfilled along all the 
trajectory, the interval dsF never vanishes and all the 
singularities in the curvatures are avoided, also for pe- 
culiar dynamical systems as the Bianchi IX models, and 
we can geometrize the full motion. The Finsler metric is 
given by 



i d 2 A 2 (r a r' a ) 
o« (x a x' a ) - - { ' ' 

where A is the homogeneous Lagrangian 



(3) 



A(x a ,x"*) = £(t,x i ,^)x'°, (a, /3, 7 = 0,..., A), (4) 

x° = t and x' a = dx a /dsF- As we see, the metric and 
all the other geometrical quantities depend on both co- 
ordinates and velocities. The mathematical tool to in- 
vestigate the stability, in the geometrical framework, is 
given by the geodesic deviation equations which, in local 
coordinates, read 



V 

ds 



Vz a 



R a 



dx p „dx s 



-WV 7 . ( 5 ) 



where z a is the perturbation to a geodesic, V/ds is the 
total (covariant) derivative along the geodesic, R a p 7 s is 
one of the curvature tensors defined in Finsler Geometry 
p7t and TL is the so-called stability tensor fl2] , B)]. 
The dynamical behaviour of the system is led by this 
tensor, whose eigenvalues, {Ac}, are the sectional curva- 
tures and whose trace Ric is the analogous of the Ricci 
curvature (along the geodesic) in Riemannian geometry. 
In analogy with the Lyapunov exponent, we define an 
instability exponent 



Si = lim lim 

s^oo z(0) — >0 



1 z(s) 
— In . . 

S 2(0) 



(6) 



which now turns to be a quantitative indicator, invari- 
ant under change of coordinates x a = (t,x l ), giving a 
measure of a possible exponential growth of the pertur- 
bation z(s), defined as the natural norm on the manifold, 
z 2 = g av z a z v , in terms of the geodesic time parameter. 
In order to grasp some qualitative informations on the 
stability of the flow, some simplified approximated ver- 
sions of the (N + 1)— equations (JsJ) can be derived, as the 
ones obtained substituting the covariant derivatives with 
ordinary ones and the rhs of eq. (|b|) with the N non trivial 
eigenvalues of TL: 



ds 2 



+ Aj2 (i) = (i = l,...,N) . 



(7) 
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When N ^> 1, a further simplification is attained sum- 
ming the N equations (Q), obtaining the single scalar 
equation: 



d z 



Ric 

lv* = 



(8) 



which has been successfully applied, for some TV-body 
systems @^|. 

In the case of the Bianchi IX dynamics we have consid- 
ered the Lagrangian 



C = K — U(a, /?+,/?_) 



(9) 



where 2K = j3i_ + /3_ — A , C/ is the minisuperspace 
potential @, 3a = ln(a6c), 6/3_ = V3 ■ ln(a/6) and 
6/3+ = ln(afe/c 2 ), a, & and c are the scale factors and 
the dot means differentiation with respect to the time 
parameter dr = dt/abc, defined in terms of the cosmo- 
logical proper time t. The time component of the field 
equations imposes the constraint H = K + U = 0, whose 
fulfillment must be checked carefully. We numerically 
integrated the equations of motions, using two different 
methods: a fourth order Runge-Kutta scheme and a sym- 
plectic algorithm, both with a variable time step. The 
quantity A = \H\/(j3 2 _ + fit + a 2 ) has been monitored in 
order to test the accurancy of the two algorithms pjj and 
we find that it never exceeds the value 10~ 5 . In fig.|l| the 
three non vanishing sectional curvatures and the Ricci 
curvature have been reported versus sf- All them evolve 
towards zero but some peaks are present in coincidence 
with the bounces on the potential walls. Note that Ai^ 
are not sign definite but both are negligible compared to 
A3, which is instead always positive. Therefore, in con- 
strast with the Jacobi metric, also the Ricci curvature 
is always positive. It has been already remarked that in 
some cases equations (|7]) and (|8|) describe correctly the 
qualitative features of dynamics, | fl9|| , and some authors 
fl3"| , [22f explicitly considered them. So, in fig.|| we plot 
the instability exponents calculated from these equations. 
Note that all the Si go to zero, confirming the absence 
of any global exponential instability. However, to be rig- 
orous, we also considered the exact geodesic deviation 
equations in the following way. As the two variables r 



and a are monotonic functions |23| of Sp, the possibly 
unstable behaviour cannot be ascribed to them, so we 
restrict ourselves to consider the growth of the projec- 
tion z of the perturbation z in the plane /?+/3_, which is 
also the physically meaninguful space in this context, 



z 2 



(10) 



These are the only variables which can give rise to an ex- 
ponential instability, as they describe the anisotropy of 
the model. The metric in this submanifold is moreover 
positive definite. In fig.^j we plot the instability expo- 
nent for z for several initial conditions. As in the case of 



fig||, for all the conditions we consider, there is a clear 
absence of any exponential divergence. However, it must 
be stressed again that, in correspondence of the bounces, 
the Si tends always to increase (see the inset in fig.||), 
though the amount of this rise reduces during the evolu- 
tion; correspondingly, the sf~ interval between two suc- 
cessive bounces increases, roughly exponentially, so that 
the instability exponents undoubtedly decrease to zero as 
Si cx Sp 1 . From the strict definition of Chaos, we should 
then conclude that the system is not chaotic due to the 
absence of a global exponential divergence of nearby tra- 
jectories. However, the local sharp increase of the Si's 
reveals that a different kind of impredictable behaviour, 
confined in very short time intervals, could be present, as 
also outlined by Contopoulos et al. |24|, though within 
a very different approach. The unpredictability in the 
duration of each long era and in the number of Kasncr 
epochs for each era, is another signature of the stochastic 
nature of the Bianchi IX dynamics, which reflects itself 
in the sharp rises, due to the chaotic reflections from the 
potential walls. It is obvious that stochasticity depends 
on the definition adopted and also on the observable cho- 
sen. 

In summary, we presented a gauge invariant approach to 
Chaos for Bianchi IX dynamics, without assuming any 
kind of approximations. Our results reveal the absence 
of any global chaotic behaviour even if they suggest the 
occurrence of a chaotic scattering phenomenon. Further 
details and a wider discussion on the subject will be pre- 
sented elsewhere. 
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FIG. 3. Evolution of the instability exponents 5 



0+B- 



de- 



fined for z in equation ([to]), for several initial conditions. In 
the inset it is plotted the curve corresponding to the same 
initial condition of figsjl| and || 



FIG. 1. Evolution of the sectional curvatures Ai, A2 and 
A3 and of the Ricci curvature per degree of freedom, Ric/3, 
for the initial conditions (3+ = 0.33, 0+ = —0.23, j3- = 0, 
$- = —0.64 and a — —0.3. Note that only sf is plotted in 
log scale. 



FIG. 2. Evolution of the instability exponents 5^1, S\ 2 , 5\ 3 
and ^Ric/s for the approximated geodesic deviation equations 
(0) and (ji|). The initial conditions are the same as in fig.[j]. 
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